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Abstract Wc consider the one dimensional symmetric simple exclusion process (SSEP) with additional births and 
CNl ' deaths restricted to a subset of configurations where there is a leftmost hole and a rightmost particle. At a fixed 
Vh ! rate birth of particles occur at the position of the leftmost hole and at the same rate, independently, the rightmost 
C*H' particle dies. We prove convergence to a hydrodynamic limit and discuss its relation with a Stefan problem. 
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r ~7 1 ' 1 Introduction 

P< 

Ph . The Stefan problem is a free boundary problem. In the simplest case the evolution is given by the linear heat 
equation in a domain D, which itself changes in time with a law which depends on the same solution, for instance 
the local velocity of a point of the boundary is proportional to the local gradient of the solution. 

The purpose of this paper is to study a particle version of such free boundary problems. The linear heat equation 
is in our case replaced by the one dimensional symmetric simple exclusion process, SSEP, in the set of configurations 
having a rightmost particle and a leftmost hole. Furthermore, at a given rate the rightmost particle dies and at the 
same rate, independently, a birth of a particle occurs at the position of the leftmost hole. The leftmost hole and 
the rightmost particle are the boundaries of the system; for this reason we call it free boundaries. 

Call r\ £ {0, 1} Z a particle configuration, think of r\ as the subset of Z occupied by particles and consider those r\ 

pC ' having a rightmost particle located at 11(77) := max(?7) and a leftmost hole located at h{n) := min(Z \ if). Let r\ t be 
the Markov process performing symmetric simple exclusion process at rate i and such that the rightmost particle 
and the leftmost hole are killed at rate J: 



V —> V \ i 1 ^ 77 )} an d V —* V U {L(7y)} at rate J each. 



Since particles arc injected to the left and extracted from the right, J can be seen as the average current of particles 
through the sistem. It is well known that under a diffusive space and time scaling the collective behavior of the 
SSEP is ruled by the linear heat equation j2j . We perform the same scaling with a parameter e such that time is 
r> ■ e~ 2 t, space e~ 1 r and the killing is J = j(e) = ej, where j is the macroscopic current. 

Consider a function p : R — > [0, 1] identically zero to the right of R,(p) := sup{r : p(r) > 0} < oo and identically 
one to the left of h(p) :— inf{r : p(r) < 1} > — oo and continuous in (l(p), r(/o)). Call 1Z the set of functions with 
those properties. We consider a macroscopic density p € 7Z and ask the initial configuration rj 1 - 6 ' indexed by e to 
approach the density p as follows: 



f b 
lim sup \e \ ip £ '(x)— / p(r)di 



(1.1) 



sx£[a,b] 

Our main result is 

Theorem 1. Let rff be the process with killing at rate je and with initial configuration -q 1 - 6 ' satisfying (jl.ip . Then 
for t > there exists a function p t <E 71 such that po = p and 

limPfsupe Yl Vie-'( x )- Pt{r)dr > 7] = for all 7 > 0. (1.2) 



We do not have a description of the evolution of pt , but we can construct it as limit of approximate evolutions 
p\ as 5 > goes to zero. The functions p\ satisfy a "discretized Stefan problem" , namely p\ evolves according to 
the heat equation in the intervals [nS, (n + 1)5) and at times n5 takes the rightmost portion j5 of mass from the 
right and puts it on the left. More precisely, define the <5-quantiles R 5 (p) and L 5 (p) by 

p{r)dr = 5, / (1 - p{r))dr = 5, (1.3) 



'R s (p) J -co 

Define also 

{1 iir<L 5 (p) 

p(r) ifL S (p)<r<R S (p) (1.4) 

ifr>R s (p), 

and let G t (r,r') be the Gaussian kernel ( see (|5.10[) ). Set iteratively 

s . = iGt-nSP S nS , i£t£[n5,{n + l)5), n = 0,l,... 
Pt " '{&"&>-> iit = n5, n = l,2,.... 

which is well defined for 5 small enough. We prove 

Theorem 2. Let p G 1Z and p t be the evolution of Theorem^ with initial datum, p. Let p t be the evolution (| 1 .5[) 
with the same initial datum. Then for any a < b real numbers and for any sequence 5 n so that 5 n +\ = 5 n /2 we have 

b ,-b 

p s t "(r)dr- p t (r)dr < 2j5 n , Vt > (1.6) 



Theorems Q] and [2] are proved at the end of Subsection 16.11 

A formal limit of the above discrete Stefan problem leads to conjecture that p t solves 

R ,L ,p(r,0) given 

p(U,t) = l, p(K t ,t)=0; ^-(j Jt ,t) = ^-{K t ,t) = -2j 

or or 

where R t := R.(pt) and L t := L(p 4 ). (11.71) is complemented by equations for the motion of the boundaries: 

dh t 1 d 2 p dK t 1 d 2 p 

^-Y 3 ^ {U ^ ^W^ (L8) 

A similar equation has been derived in [7] for a different process, see also reference therein for particle models for 
the Stefan problem and [S] for the derivation of a stationary Stefan problem when phase transitions are present. 

The model we study in this paper is inspired by previous works on Fourier law with current reservoirs, [3]-[4]; 
its actual formulation came out from discussions with Stefano Olla to whom we are indebted. 

In Section [2] we define the particle process and in Section [3] we prove the existence of a unique invariant measure 
for the process as seen from the median. 

Due to the non local nature of the birth-death process the usual techniques for hydrodynamic limit fail. To 
overcome this problem we use inequalities based on imbedding the particle process in an interface process. The 
relationship between the one dimensional nearest neighbors simple exclusion process and the interface process is 
known since the seminal paper by Rost [9j , where he established the hydrodynamics of the the asymmetric simple 
exclusion process. In Section |4]wc define the interface dynamics and show that our particle process can be realized 
in terms of the interface dynamics. We also introduce the delta interface processes which correspond to the delta 
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L M R 

Figure 1: A typical configuration in X, black and white circles represent respectively particles and holes. R is the 
position of the rightmost particle, L the position of the leftmost hole and M is the median. 

particle processes where the killing of particles and holes are grouped together and occur only a finite number of 
times (uniformly in the hydrodynamic limit) and give as limit the delta macroscopic evolution defined in (If .51) . 

In Section @] we also establish basic inequalities between the true and the delta interface dynamics by giving a 
simultaneous explicit graphical construction of all of them; that is, a coupling. In Scction[5]wc prove convergence in 
the hydrodynamic limit. The proof uses that any limit point of the true interface dynamics is squeezed in between 
two approximate evolutions which depend on an approximating parameter 5 > 0. In Section [5] we establish basic 
properties of the macroscopic evolution. In particular, we prove the existence of a stationary solution for the limit 
evolution and use this result to prove that at any positive time the particle density (in the limit evolution) is 
identically and I outside of a compact. In the last section we summarize the results. 



2 The free boundary SSEP 

The space of particle configurations is 

X := { V £ {0, 1} Z : £ r,{x) < oo, £(1 - r)(x)) < oo} 

x>0 x<0 

that is, for any r\ £ X the number of particles to the right of the origin and the number of holes to its left are both 
finite. A configuration r\ £ X has a rightmost particle located at r(^) and a leftmost hole located at h(rj), where 

r(ji) := max{x £ Z : t)(x) = I}; l(t/) := min{a; <G Z : r](x) = 0}. (2.1) 

We define the median M(rf) of a configuration r\ £ X as the unique meZ+| such that 

5>(z) -£ (1 -??(*)) =0, (2.2) 

x>m x<m 

that is, the number of particles to the right of M(r)) is the same as the number of holes to its left. The definition is 
well-posed because for r] £ X, as m increases by one, (|2.2|) increases by one and goes to ±oo as m goes to =Foo. We 
next define a family of dynamics indexed by a current J > 0. The particle dynamics is a (countable state) Markov 
process on X whose generator is 

Ll tt =L +L ] R + L[, (2.3) 

where 

WW^E^/f^V/W] (2-4) 

with r) x ' x+1 (y) = rj(y) if y j^ x,x + 1, rf' x+l {x) = r)(x + I), r] x ' x+1 (x + I) = 77(21); and 

L[f(v) := j(/faU{Lfa)}) - /(??)); L J J( V ) := j(/(rA{R(f/)}) - /(»?)), (2-5) 
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L M X R 

Figure 2: The effect of three SSEP jumps: the upper line is before and the bottom line is after the jumps. The 
jumps are between the positions L and L + 1, M — 5 an( l M + 2 an d x and x + 1. None of these jumps change the 
position of M. 
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Figure 3: The effect of killing the rightmost particle: the upperline is before the killing and the bottom line is after 
the killing. The position R of the rightmost particle moves to the left by 3 (for this configuration) and M moves by 
1 to the left (for any configuration). Analogously, the killing of the leftmost hole moves M one unit to the right. 

where n is identified with the set of occupied sites {x € Z : i](x) = 1}. 

In other words, particles perform symmetric simple exclusion with generator Lq and, at rate J, the rightmost 
particle is "killed" and replaced by a hole and at the same rate independently the leftmost hole is killed and a 
particle is born at its place. We omit the proof that the process is well defined at all times. Denote by rjt the 
process with generator L? t . 

Denote by B t and A t the number of particles killed, respectively born, in the time interval [0,i\. By definition, A t 
and B t are independent Poisson processes with rate J. 

Lemma 3. For any r)o € X, 

M( Vt ) = m(? ? o) + A t - B t (2.6) 

That is, the marginal distribution of the median of r\ t is a continuous time symmetric nearest neighbor random walk 
on Z + i with rate J to jump right and J to jump left. 

Proof. Note that (a) jumps due to the exclusion dynamics do not change m(-), (b) when the rightmost particle dies, 
m(-) decreases by 1 and (c) when the leftmost hole dies (and a particle appears in its place), m(-) increases by 1. 
This shows (12^1). □ 



3 The process as seen from the median 

Since M(rjt) is a symmetric simple random walk the law of rjt cannot be tight but we show below that the process 
as seen from the median has a unique invariant measure. Let 

Vt '•= M (j7 t )-l/2??i! 

where for y G Z, the translation y : X — > X is the map (8 y i])(x) = rj(x — y). Clearly f] t E X Q := {77 <G X : M(n) = 
1/2} and we have that r) t is a Markov process on X° with generator 

LU = L ° + l i + Z 'h (3- 1 ) 

where for any i] G X a , 

L\f{n) := j(/(ft_ x fa \ {K(r,)})) - /fa)); L[f(n) := j(/(^ U { L fa)})) - f( V )}) 



Before stating the result we introduce some notation. For 77 £ X, define 

N°( V ):= E 0--v(x)), N\ V ):= ]T r,(x) (3.2) 

x<r(ti) x>L(r/) 

the number of holes in 77 to the left of the rightmost particle and the number of particles in 77 to the right of leftmost 
hole, respectively. Clearly 

N° + N 1 = r-l + 1. (3.3) 

Let 77° £ X° be the Heaviside configuration given by 

77°(a;) = l{x < 0}. (3.4) 

Define \j) on X by 

1>(jj) := J2 (1 - V(?))V(V) (3-5) 

x<y 

Observe that ip(r)) < 00 for all 77 £ X. If 77 £ X°, ip(n) is the number of jumps 10 — > 01 needed to get from 77 to 77. 
In particular, ip(r] ) = 0. 

Theorem 4. For any J > the process fjt has a unique invariant measure fji 3 on X° and 

^[r-L + 1] = 1 (3.6) 

Proof. We show that ip is a Lyapunov function for fj t by computing L\ vt ip. Since ip(9 x rj) = ip(r]) for all x and 77, 

(Z£ + LiMr,) = (L[ + L^{rj) = -j(N°( V ) N^rj)) = -.1(11(77) - Lfa) + l) (3.7) 

the second identity holds because when the rightmost particle disappears, ip decreases by A^ , the number of holes 
to the left of R. Analogously, when the leftmost hole disappears, ip decreases by TV 1 , the number of particles to the 
right of L. The third identity follows from (|3.3j) . 

Since ip increases by one when there is a transition 10 — > 01 while it decreases by one due to the opposite 
transition: 

Lotto) = \ X>(*)(1 - V(x + 1)) - (1 - V(x))r,(x + 1)} = i, (3.8) 

X 

because for any configuration 77 £ X, the number of pairs 10 exceeds by one the number of pairs 01. The sums in 
are finite for 77 £ X. 
Call v t the law of 77^ starting from a configuration 77 £ X°. For any t > we have v t ip < 00 and 



1 , rt 

■{v t 1p - V Q lh) = - I 

1 ft 



1 1 r 1 j r 

(v t il> - iwl>) = - / v s [U^}ds = --- I Vs [b. - l + 1] ds. (3.9) 

Wo *• t Jo 



Hence, calling fi t := \ J v s ds, 

M t [R-L + l]<— + -^- (3.10) 

Since for any c we have #{77 £ X° :R — L + l<c}< 00, then the family {fj,t,t > 0} of probabilities on X° is 
tight on X° and it has therefore a limit point fi. fi is stationary by construction and unique because the process is 
irreducible. 

Identity (|3.6|) follows from (|3.9p by replacing v t by /j, but we need first to show that flip < 00. Introduce a new 
Lyapunov function ipi on X by setting 



Mv)-= J2 [l-r } (x)][l-r,{y)]r,{z) 



x<y<z 



Then, 



l V> 2 (t?) - \ ( E I 1 ~ *7(*)Mv) [i - ? ?(y + 1)] - E t 1 - ^)] t 1 - »?(jOMv + 1; 

x<y x<y 



X<J/ 



and 



l 3 Mv) = -j E [i -*/(*)][! -»?(y)] = -^°W[^°W-i] 



2 

x<y<R 

l[Mv) = -j E ^-^GOM*) = - j E[ 1_7 7( y )M z ) +J E ^) 

= -Sipirfj+JN 1 ^) 
Let ^* as defined before (|3.9[) and /it after it. Then ^^2 < o° an d since J^ [1 — r](x)]rj(x + 1) < N°(r]), 



1 / , . \ J /■' 



iV° , „i A^ 



ds 



/ \ i / r iV i v 

t [Ufa - v^2) <--l v s [- — +^ - N 1 + —{N° -1) 

so that, since N°(N° -l)> 0, 

1 riV ,i 

ptip < — t'oV'2 + Mt "t; 1" -<* ^ Constant, 

by (f373|) and (|3T0| . Then pip < oo and setting z^ = /x in ([3^9]) we get (|3~6| . D 



The theorem says that under the invariant measure the average distance between the rightmost particle and the 
leftmost hole is (2j) _1 . This is in agreement with the Fick's law (the analogue of the Fourier law for mass densities). 
In fact Fick's law states that the stationary current J flowing in a system of length £ when at the endpoints the 
densities are p± is: 

T _ 1 P+ - P- 
J 2 £ 

1/2 being the particle mobility. In our case J = 2j, p+ = 1 and p_ = hence £ = (2j) _1 . The validity of Fick's 
law in our case is however not completely obvious as the endpoints 11(17*) an d l (t7*) depend on time. 

4 The interface process 

The interfaces. Let the set of vertices be 

V := {v = (vi,v 2 ) e Z x Z + with i>i + u 2 even}. (4.1) 

and for a vertex w = (v\, v-z) € V let V v : Z — >■ Z be the cone with vertex w defined by 

K(a;):=|a;-ui|+V2. (4.2) 

Define the space of interfaces as 

3>„ := {(eZ z : \£{x)-Z(x + l)\ = l, xeZ; #{x : £{x) ^ V v (x)} < oc} 
y := U uev ^ (4.3) 



That is, an interface in y v coincides with the cone V v for all but a finite number of sites. Interfaces share the 
property "x + £(x) even" which is conserved by the dynamics defined later. For an interface £ G y define 

L (£) := sup{x e Z : £(x - y) = £(x) + y for all y > 0} (4.4) 

r(£) := inf{a; G Z : £(x + y) = £(x) + y for all y > 0}. (4.5) 

Any interface £ coincides with a cone V v outside the finite interval (l(£), r(£)). The vertex v = v(£) is the following 
function of R = R(£)), l = l(£), £(l) and £(r): 

/as £( L )^£( R ) l + r _. £(l) + £(r) l-r ,, . 

"i = "i(fl:= - 2 +~2~' «2=« 3 (0-= - 2 +^~- ( 4 - 6 ) 

Correspondence between interfaces and particle configurations. Given an interface £ we say that there is a particle 
at x if £(ie + 1) < £(x) and that there is a hole at x if £(x + 1) > £(x). This defines the map D : y — > X given by 

^) S fl(0(x) = i- ^ + 1 ^-^ ) (4.7) 

The map is clearly surjective but it is not injective as D is invariant under uniform vertical shifts: !?(£ + n) = -D(£) 
for all n £ Z. However the extremes of the non conic part of £ correspond to the leftmost hole and the rightmost 
particle of D(£) and the absise of the vertex of the cone containing £ corresponds to the median of -D(£). More 
precisely: 

Lemma 5. For any £ <G y, 

l(Z»(0)=l(0, R(D(0) = R(0 - 1. M(D(0)=« 1 (0"l/2. 

Proof. The first and second identities follows directly from the definitions. The third identity is trivially satisfied 
by any cone and its mapped particle configuration: Vi(V/ a m) = a and m(-D(V( &))) = a — 5. Any interface £ in the 
cone V( ai f,) can be attained from V( 0] &) by a finite number of moves of the type: 

£ = (...,z,z-l,z,...)->{...,z,z + l,z,...) = ? (4.8) 

The corresponding particle moves are 

r, = (.. ., 1,0,. ..)->(■•-, 0,1,. ..) = // (4.9) 

and 

D(£) = 77 if and only if D(£') = ?/. (4.10) 

Observing that u(£) = i>(£') and that M(D(£)) = m(£)(£')), we conclude that £>(£) has median a— h — i>i(0 for all 

Interface dynamics. Let £t be the Markov process on y with generator 

£L« := ^a + 4 + 4, with 
£a/(£) := ^ £{/(£ + A x /(£)}, (4.11) 

4/(0 := j[/(max{£,K (s)+( _ 1 , 1) }) - /(£)], £[/(£) := j[/(max{£, K( ? )+(i4)}) " /(0]> 

where A x £(y) := (£(x + 1) + £(x — 1) — 2£(x)) l{y = x}. The jumps of £(x) due to £a occur only when £(x) has the 
two neighbors at equal height. The jump is up by 2 if the neighboring heights are both above £(x) and down by 2 
if they are both below £(x). C 3 R acts by changing the rightmost downward variation of £ into an upward one with 
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Figure 4: The thick red line represents the interface £t; its corresponding particle configuration is r\t = -D(£t)- The 
narrow blue line represents £o, the interface at time 0; 770 = D(£o)- Two particles and one hole have been killed 
and several particles have moved due to the exclusion dynamics; in particular the third particle moved to the place 
originally occupied by the second one and the second hole has moved three units to the left. Due to the killings the 
vertex of the cone containing £ t has moved by 2(— 1, 1) + 1(1, 1) which amounts to three units up and one unit left 
so that Co e 3^(0,0) while £ t <G ^(-1,3)- 

the interface to its right being a straight line with slope 1. The cone containing the updated interface is obtained 
from the previous cone by a translation up by 1 and left by 1; see Figure [U A symmetrical picture describes the 
action of £ J L . We denote by 



M = #0 um P s due to £ J L in [0, £]); B t = # (jumps due to C 3 R in [0, t]) 

A t and B t are independent Poisson processes of intensity J. 

The interface evolution induces via the map D the particle evolution described in Section [2] 



(4.12) 



Lemma 6. If the interface process £t has generator /Cf ntor , then the particle process r\t '■— -D(£t) is Markov with 
generator Z/ art , defined in (|2.3|) . Moreover, if £0 S 3^(o,o); then 



Ct(0) = 2B t + 2^2r, t (x) 



(4.13) 



x>a 



Proof. It just follows from the definitions of £t and D that D(^ t ) is Markov with generator L J art . Calling R = 
R (6) = Hvt) + 1 we get from (|4~7|) 

x>0 

From ([46]) . ^(r) = R - (vi - v 2 ) and since V\ = A t - B t and v 2 = A t + B t , we get (I4.13J) . D 



Harris graphical construction Wc construct explicitly the interface process £( as a function of the initial 
interface and of the Poisson processes governing the different jumps. 

We construct first the process with generator £a and later use it to define the process with the moving boundary 
conditions. The probability space (f2, P) is the product of independent rate-^ Poisson processes on K + indexed by 
Z x {f, 1}. A typical element of il is u = (w|, wj, x G Z). The Poisson points in w| and respectively u4 are called 
up-arrows and down-arrows, respectively. 

Wc define operators T t : f2 x ^ — > J 7 with t > where T t (uj,£) is the process with initial interface £ using the 
arrows of uj, as follows. We drop the dependence on u and write just T t £, instead. We can take ui such that at most 
one arrow occurs at any given time. 

Set To£ = £. Assume that £ t / := T t /£ is defined for all t' e [0, s]. Let t be the first arrow after s belonging to 
w| U wj: for some x such that £ s (a; + 1) = £, s (x — 1). Since £ € J?, there are a finite number of such x and t — s > 
a.s. These are the arrows involved in the evolution at time s. 

Set 

7K = T s £ for *' G [a, t) 

and (1) If t is an up-arrow, then the interface at x is set to £t(x — 1) + 1 no matter its value at t— and does not 
change at the other sites: 

T£(x) := T t -Z(x - 1) + 1; I^fo) := T t _£(y), y ^ x. 

(2) Analogously, if i is a down-arrow, 

T£{x) := T t _^x - 1) - 1; T t £(y) := T t _£(y), 2/ ^ s. 

The reader can show that the process T t £ so defined is Markov and evolves with the generator Ca with initial 
interface £ at time t = 0. 

In the next definition we need to use the operator T t in different time intervals. With this in mind wc define 

T M (w,S) :=Tt- a (8- a w,£), (4.14) 

where 8 s uj is the translation by s of the arrows in u. That is, TL t i has the same distribution as T 4 _ s but uses the 
arrows in u belonging to the interval [s,t\. We drop the dependence on uj in the notation and write simply Tr st i£. 

Generalizing the boundary conditions Consider the partial order in the vertex space V given by 

v < v' if V v (x) < V v >(x) for all i€Z. (4.15) 

So that vertex order corresponds to cone order. Let z t G V be a non decreasing path of vertices with finite number 
of finite jumps in finite time intervals: \\z t — Zf\\ < oo for all t < t' . Let 3T s ,t] be the family of random operators 
governing the £a motion, defined in (|4.14[) . We assume that the process z t is independent of the random operators 
T[ B t] ■ Typically z t will be a function of the Poisson processes A t and B t of rate J which are independent of the 
arrows. We abuse notation and call P the probability associated to u>, A and B. 

Let = sq < si < . . . be the times of jumps of z t . Define iteratively T££ := £ for t < and 

fmax{T, z £,V Z A i£t = s n , , 

T t z £,:= I X '"^' * ; ' (4.16) 

[ T [s n ,t] T a n € lf l £ ( s n, S„+l). 

So £t evolves with £a in the intervals (s„, s n +x) and at times s n is updated to the maximun between the interface 
at time s n — and the cone with vertex z Sn . 

Given an initial interface £ € y and a number 5 > 0, wc will consider the following choices for Zt, denoted by 
Ot, Rt and z t ~ and z t : 

0/ 

i?/ 



«,- 



w(0, for all t; (4.17) 

w^J + ^t-Bt.At + Bt) (4.18) 

i?„5, ^' + := i?(„+i)a, te[ni,(n+l)i),n>0. (4.19) 



When the path is Ot, the cone does not move and the resulting process has generator £&■ When the path is Rt the 
process has generator £j ntcr . The process with path z t '~ records the increments of Rt in the intervals [n<5, (n + 1)5) 
and takes the maximun of the interface at the end of this interval and the cone with vertex R n $ . The process with 
path z t records the increments of Rt in the intervals [nS, (n + 1)8) but takes the maximun between the interface 
at the beginning of this interval and the cone with center Rt n +us- The processes with paths (|4. 10[) will be used 
later. 

Monotonicity and attractivity Consider the natural partial order in y given by: £ < £' if and only if £(x) < 
£'(x) for all x <E Z. Use this order to define stochastic order for random interfaces in y. 

Let £j and £ 2 be two realizations of a stochastic process £ t on y with initial random interfaces £ x and £ 2 , 
respectively. We say that the process £< is attractive if the following holds: 

If £} < £ 2 stochastically, then $ < £ 2 stochastically, for all t > 0. (4.20) 

Lemma 7. If £ < £' then T t £ < T t £' almost surely. As a consequence, the process with generator Ca is attractive. 

Proof. Consider £ < £' and call & = T t £, $ = T t £', assume that &_ < £_, &_(x± 1) = £ t '_(x± 1) and that t e w|. 
Then £*(x) = £' t (x) = £t-(x — 1) + 1 no matter the values £t-(x) and ££_(x). Analogous argument applies when 
t E u^.. Hence an arrow does not change the order if the two interfaces coincide at x ± 1 . If at least one of the 
neighbors of x in £ t _ is different of the corresponding neighbor in £ t '_, a similar argument shows that no jump can 
break the domination. We have proven that if the the £ process is dominated by the £' process just before an arrow, 
then the domination persist after the jump(s) produced by the arrow. Since the set of involved arrows is finite in 
any time interval, an iterative argument concludes the proof. □ 

Remark It is usual to realize the process with generator £a by introducing only one rate- ^ Poisson process of 
marks uj x ,x € Z associated to each x and updating £(x) — > £(x — 1) ± 1 7^ £(x) whenever a w^-mark appears 
provided £(x — 1) = £(x + 1). This is indeed a realization of the process, but order is not preserved: the jumps 
cross if two interfaces coincide at x ± 1 and differ by 2 at x at the updating time of x. Using the up- arrows and 
down-arrows, only one of the interfaces jumps, and order is preserved. 

Lemma 8. If £ € y, then v(T t t;) = v(£) a.s. As a consequence, V v r^ < T t £. Furthermore, for any non decreasing 
path z t such that i>(£) = zq, u(T t 2 £) = z t , which implies V Zt < T" t z £. 

Proof. By the definition of T t £, any time L jumps due to an arrow, the opposite jump is performed by £(l) and 
analogously, any jump of R is replicated by £(r). This implies 

L(T t + tfL(T t 0) = L(£) + f(L(0), R(TtO - f (R(T t 0) = R(0 - £(R(0) ( 4 -21) 

Putting these identities in the definition (|4.6[) of v(£) and noting that the total number of jumps of L and R is a.s. 
finite (it is dominated by a Poisson process of rate 2), we get v(T t £) = i>(£). 

The fact that v(T t 2 £) = z t is true by definition at the z-events and by the first part of the lemma, it is true for 
t G [s, s') where s and s' are successive z-events. D 

Proposition 9. Let Zt and z' t be non decreasing paths in V. Then, 

if z t < z' t for allt>0 and £, < £', then T 2 £ < T t z ' £ for all t > a.s.. (4.22) 

In particular, by taking z = z', Proposition [9] says that for any non decreasing path z t , the process T/£ is 
attractive. 

Proof. Since the operator Tr a t i is monotone by Lemma [71 we only need to check that the inequality T t z £ < T t 2 £ is 
preserved when z and z' events occur. We thus suppose the inequality is satisfied for all s < t and this is evidently 
still true if t is a z'-event. If instead t is a z-event, T*_£ < T*_£' implies 

T t z i = max{T^,U Zt } < max{T t W , V Zt } < max{T t W, V z[ ] = T? £, 
where the last inequality follows from V Zt < V z > and the last identity from the last inequality in Lemma [5J □ 
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The interface process as seen from the vertex Take o = (0, 0) and £ £ y o and define the interface process 
Tft; as seen from the vertex by 

Tf£(x) = 6 at T?Z(x), (4.23) 

where the translation by v = (i>i, 1*2) G V on y is defined by v ^(x) := £(x — wi) — U2. Of course T t 2 £ G 3^ . Also, if 
a vertex u g V satisfies u > o (that is, V v > V ) and £ G 3^ OJ then 

0„£ < £ and max{6W, V } < f . (4.24) 

Monotonicity We show that if the initial interface £' dominates £ and any jump of z' is replicated by z (that 
is, z "jumps more" than z'), then the interface process as seen from z' dominates the one as seen from z. More 
precisely, 

Proposition 10. Let z t and z' t be non decreasing paths on V. Then, 

if z t - z t - > z' t - z' t _ for allt>0 and £ < £', then fft, < f/'f for all t > a.s.. (4.25) 

Proof. Since by Lemma UJ the domination is preserved in intervals with no z or z' events, it suffices to take care of 
those events. Assume that T t z _£ < T t z _£' and that t is an z' event, then t is also an z event by hypothesis and 

f# = max {6 Zt - Zt _f t U, V ) < max {6 Zt - Zt _f t W, V ) = 7?'? 

If t is an z event and not an z' event, then 

ff£ = m a x{9 Zt - Zt _TU,V } < ff_t < f t W < f t W = ???. (4.26) 

where the first inequality in (|4.26[) follows from (|4.24[) . D 

Let Zt(l) — Zq(j) := Rt = (A t — B t , A t + B t ), where A t and B t are independent Poisson processes of rate J. For 
each J > the interface process £ t = T t £ has generator is £J ntor given by 

4 tor = £ A + 4+4 (4.27) 

where £a was defined in (|4.11j) and the other generators govern the updating with the maximun of the cone and 
the corresponding translation of the origin: 

4/(0 ^[/(maxl^^K}) -/(£)], 
£1/(0 := j[/(max{0 (ljl) £, V }) - /(£)] . 

For £ <G y o , the process D(T*£) has the same law as the particle process f\ t defined in Section [3] with initial particle 
configuration D(£) € Xq. The map D : y^ — > Xq (interfaces with vertex in the origin to particle configurations 
with median — ^) is bijective. Since the process fj has a unique invariant measure on Xq, l~* has a unique invariant 
measure jx } on y$ . Taking z — z' in (|4.25|) we get immediately the following. 

Corollary 11. The process £t := T t £ is attractive. Furthermore, if £0 = Vo the law of ^t is non decreasing in t, 
is stochastically dominated by the invariant measure fi 3 and converges to ft? as t — > 00. 

Corollary 12. The invariant measures fj, 1 for the interface processes T t £ are stochastically ordered: 

If J > J'j tten /5 J < /2 J 
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5 Hydrodynamic limit 

It is well known that in the diffusive scaling limit (space scaled as e _1 and time as e -2 ) the SSEP process alone 
converges in the hydrodynamic limit to the linear heat equation: 

dpt ld 2 p t . , 

¥ = 2^' P0 = p (5 - 1} 

Since the distance between the first hole and last particle is random, it is not clear a-priori that our model should be 
scaled diffusively as well. From Theorem [4] we know that at equilibrium the mean distance between the rightmost 
particle and the leftmost hole is of order J" 1 and this suggest that together with the above diffusive scaling we 
should scale J proportionally to e. Indeed we will prove that under such scaling limit the density of the process 
converges as e — > to a deterministic evolution. 

5.1 Results 

Hydrodynamics of interfaces. Initial configurations. For each e > the interface evolution starts from an 
interface ^ s ' such that: 

lim sup H {6] (x) - <j> {ex)\ = (5.2) 



E->0 



j-e- 



where <t>o( r ) = r f° r au r > R (<W > 0; 0o( r ) = — r for all r < h(<p) < and cf>o is differcntiable in (l((/>),r(</>)) with 
derivative 4>b such that sup |0o( r )l < 1- 

re(L(0),R(0)) 

Fix the macroscopic current j > 0, define j(e) := je and call 

Q e := interface process with generator Ca + £l £ + £[ £ , starting from ^ e ' . (5-3) 

Under £ t the particles and holes are killed at rate ej. The following result is the hydrodynamic limit for the 
interface process. 

Theorem 13. There is a function 4>t(f), t > 0, r £ R, so that for any 7 > and t > 0: 



lim P 

e-H) 



suv\e£} 2t {x)-4>t{ex)\> 1 =0 (5.4) 



The proof will be given in subsection 15.51 and properties of 4>t{f) will be discussed in Section [B] 



Hydrodynamics of particles. As we shall prove at the end of subsection 16. II The hydrodynamic limit for the 
particle process Theorem Q] is a corollary of Theorem [13] 



The proof that (|5.1[) is the hydrodynamic limit for the SSEP is quite simple because the correlation functions 
obey closed equations. The addition of birth-death processes spoils such a property but if the rates are cylinder 
functions and are "small" (births and deaths happen at rate e 2 ) the proofs carry over and the limit is a reaction 
diffusion equation as in [T] . In our case the birth-death rates are not as small (because the killings rate are of order 
e) but the main difficulty is that the killings are highly non local functions of the configuration, since births and 
deaths occur at the position of the leftmost hole and the rightmost particle, respectively. This spoils completely an 
analysis based on the study of the hierarchy of the correlation functions. 

The way out is to use inequalities namely to sandwich the interface process between two delta processes, using 
Proposition 1101 The corresponding particle delta process behave as the exclusion dynamics in macroscopic time 
intervals of lenght S and the (accumulated) killings occur at the extremes of those intervals. 
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5.2 The delta processes 

The delta interface process £ t ' is defined via (|4.16[) and (|4.19j) by 

&* := T t z± £, with z± := zf & (5.5) 

This process is obtained by patching together finitely many pieces of the evolution with generator £a as explained 
after gUJ). 

The delta particle processes are defined by using the map (|4.7[) . setting r\ t = D(^ t ). In other words, r) t '~ 
evolves with the generator Lq in the intervals [nS, (n + 1)5), and rj^ is obtained from rj n s- by removing its B$ 
rightmost particles and its A$ leftmost holes, where A t and B t are the independent Poisson processes with intensity 
J. The interpretation of r) t is analogous but the removal of particles and holes is done at the beginning rather 
than at the end of each time interval. For a particle configuration ry € X and positive integers a and b define the 
quantiles L a (?y) and R b (rj) as the lattice points satisfying 

J2 v{x) = b, Y, 0-~v(x))=a. (5.6) 

x>n b (ri) x<L a (r/) 

and define 

r a ' b ( V )(x) = V (x)(l - w) + (1 - v(x))l x <^; (5.7) 

this is the configuration obtained from r\ by erasing particles and holes as explained above. 

We abuse notation and write Ty st \r\ the evolution D(TL t i£), with 77 = D£. Then the delta particle processes 
satisfy the following 

Vt 

Vt' (5-8) 

<5.+ 

Vt 
where recall At and Bt are independent Poisson processes of parameter J. 

The rescaled delta processes. We consider a family of processes indexed by e by considering s~ 2 5 instead of 5 
and ej instead of J. We call £*' ' and 77*' ' the interface and particle processes so obtained. Later we consider 
those processes with time rescaled by e~ 2 t and space by e~ 1 r. 

For any fixed (macroscopic) 5 > 0, we will prove the existence of <j> t ' the limit as e — > of e£f'_' a ■ We also 
prove that <j) t ' are close to each other and that their difference vanishes as 5 — > 0. Taking 5 to zero, their common 
limit <j) t is then the hydrodynamic limit of the rescaled evolution Q e ' -as this is squeezed between £|' '~ and £f' ' + . 

While the above outline involves the interface process alone (which then implies convergence for the particle 
process as well), yet the analysis of the limit as e — >• of £j ' ' is more conveniently studied by looking at the delta 
particle process r] t ' ' and then translating the results to the delta interface process. We start from the particle 
model defining first the corresponding approximate macroscopic density delta-evolutions and then prove existence 
of the hydrodynamic limit for the delta particle process. 



T[m 


s.± 

5,t]Vn8 . 




for t G [ 


n5, (n 


+ 1)8) 


T A, 


"*"»& 


- ) 


fort 


= n5 




r A ( 


n + l)S,B( n 


.+ 1)1 


<5,+ 
Vt , 


for t 


= n5 



In this subsection we fix 5 > and define the macroscopic delta evolutions of densities p t ' and interfaces cf> t ' . 



5.3 The macroscopic delta-evolutions. 

this subsection we fix 5 > and define the macros 
Preliminary results. For any density p : R — > [0, 1] define the R+ U {+00} valued functions 

/•oo pr 

F(r;p):= p(r')dr' ', F(r;p):= (1 - p(r'))dr' 



representing the mass of p to the right of r and the antimass to its left. These are the macroscopic analogues of 
the number of particles, respectively holes, to the right, respectively left, of r. We introduce two disjoint subsets of 
densities called 1Z and U. 
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Let TZ be the set of densities p £ L°°(R, [0, 1]) satisfying the following conditions: 

(i) F(0, p) = F(0, p) < oo, that is mass to the right and antimass to the left are finite and the origin is the median 
of p. 

(ii) p has finite boundaries. That is, there exist — oo < L(p) < R(p) < oo such that p has finite support [l(p), R(p)]: 
p put no mass to the right of R(p) and no antimass to the left of L(p). 

(iii) p is continuous in the interior of the support. That is, if h(p) < R(p) then p(r) is continuous in (h(p),R(p)) 
with values in (0, 1). 

We call TZ* the set of p £ TZ such that R.(p) > 0. 

Let W = {«€ C(K, (0, 1)) : F(0, u) = F(Q, u) < oo}. This is the set of continuous densities with median 0. 
Let h £ TZ be the Heaviside density defined by h(r) = lr r<0 }. Clearly h(h) = r(K) = 0. 

Lemma 14. F(r; u) < oo and F(r; u) < oo for all r £ M. and u £ U U TZ. 

If p £ TZ then R(p) > and h(p) < 0; the two inequalities are strict unless p = h. 

IfR(p) > then h(p) < and the derivatives F'(r; p) and F'(r; p) exist in (l(/o), R(p)) where they are respectively 
strictly negative and positive. 

IfueU then F, F G C X (R) and F'(r; u) < and F'(r; u)>0 for all r e R. 

For any u E U UTZ and for any 5 > there are unique points R (u), L (u) such that 

F(l s ;u)=6; F(r s ;u)=S. (5.9) 

IfF(0;u) | <5 then R S (u) | and l/(u) | 0. 

Proof. Let u EU UTZ then for any ret 

/>oo />r 

F(r;u) = / u{r')dr — f u(r')dr < F(0; w) + |r| < oo, 
Jo Jo 

with an analogous argument showing that also F(r;u) < oo. If p £ TZ then R\(p) > because if R(p) < then 
F(0;p) = 0. Since i^ + (0;p) = F(0;p), then F(0;p) = and this gives a contradiction since R(p) < and p = 1 is 
not allowed. Moreover R(p) = if and only if p is the Heaviside density because F(0; p) = F(0; p). 

If u e U then F'(r;u) = -u(r) < and F'(r;u) = 1 - u(r) > 0, by the definition of U. If p £ TZ and 
R(p) > then h(p) < and by the definition of TZ, p is continuous in (l(/?),r(/o)) and away from and 1. Hence 
F'(r; p) = -p(r) < and F'(r; p) = 1 - p(r) > for all r e (l(/j), R(p)). 

Let u £ UUTZ. By the monotonicity of F(r; u) if F(0; u) > S then R s (u) > 0, while if F(0; u) < <5 then R s (u) < 
with the analogous property for L s (u). □ 

Definition of T s and G t . We call T 5 : U U ft ->• ^ the following map. If R s (u) < and therefore L s (u) > 
we set T s (u) = h, the Heaviside density. If instead R s (u) > and hence L s (u) < we set p = T s (u) equal to 
for r > R s (u), equal to 1 for r < i/(u) and equal to u elsewhere. In this latter case p = r ,5 (w) £ TZ* (that is, 
R(p) = R S (u) > 0). Thus V s acts by removing a portion S of mass from the right of R S (u) and put it back to the 
left of i/(u). 

Denote by G t the Gaussian kernel: 

G t {r,r'):^^=e- (r - r ^l 2t (5.10) 

And write Gtp(r) = J dr'Gt(r, r')p(r'). Recall that Gtp is the solution of the heat equation (|5.1[) with initial data 
P- 
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Lemma 15. Let p G TZUU. Then Gtp G hi for any t > 0. Moreover, calling 

U s = {ueU : L s (u) < < R S (u)} (5.11) 

for any j > there is S(j) > so that for any u E 1ZUU and 6 < S(j), Gsu G IAjs and therefore Y 3 (Ggu) G TV . 
Proof. Since pGKUW, we have ^(0; p) < oo. Then for any t > 0: 

F(0;G 4 p)= / dr / dr'G t {r,r')p{r') 

JO J-oo 

/•oo z'+oo /*oo /*0 

</ dr / dr , G t (r,r , )p(r / )+ dr dr'G t (r,r') < F(0; p) + c, 

JO JO JO J-oo 

where we used Fubini and then that J Gt(r, r')dr — 1 to bound the first term. Since Gtp G [0, 1], F(r; Gtp) < oo 
for all r (see the beginning of the proof of Lcmma Hl| . An analogous argument shows that F(r; Gtp) < oo for all r 
and t > and, being the solution of the heat equation, Gtp G C°° for all t > 0. To prove that Gtp € U it remains 
to show that F(0; G t p) = F(0; G t p) for all i > 0. For any t > we have 






(/ drG t p{r)-J dr[l - G t p{r)]) = 



1 o / \ 

because by (15. ip both terms are equal to — I _ n . The identity F(0;p t ) — F(0;p t ) = then follows from its 

2 Qr lr-u 
validity at time and by continuity. 

The last statement in Lemma [T5] follows from the following inequalities 

/•OO 

drG s u(r) > / drG s h{r) >CVS, C > 

which is larger than jd for i5 small enough. The first inequality follows from the fact that u is stochastically larger 
than the Heavisidc density h, in the sense that u is obtained from h by moving mass to the right. The last inequality 
follows from direct computations. □ 

Density delta evolutions We are finally ready to define the density delta evolutions p t ' . Restrict to S < S(j) 
as defined in Lemma 1151 take the initial density p G TZUU and define iteratively p ' ~ p and 

*,- . . \ G t-nSP S n s, t€[n6,(n + l)6) 
Pt \T^p s t C, t = nS. [ ' ' 

The evolution p t ' is defined as p t '~ but with initial datum p ' := T* (p). The evolution p t '~ is the one called p t 
in the Introduction. 

Interface delta evolutions The interface delta evolutions are defined as follows. Fix an initial interface <j> 
belonging to the cone with vertex at the origin and for t > define iteratively </> ' _ = 4>, cf> ' = max{^, V(o,8j)} an d 
for n > 0, 

<j>t '■= Gt-nS^n^, if t G [n5, (n + 1)5) 

max{(j) t CV(o tn sj)}, if t = n5. (5.13) 

4>t' + '■= max{<l>t-> v (.o,(n+i)6j)}, ir t = n5. 
We leave the proof of the following Lemma to the reader. It relates both definitions. 
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J,- r J.S,- 

h ■ 



Lemma 16. The delta density evolutions p t defined in (|5.12[) and the delta interface evolutions (f> t ' defined in 
(|5.13|) are related by 

^ ± (r)-^ ± (r')=2(r-r')-2 f pf'W" . (5.14) 



The initial data are related by 4>(0) = J_ (1 — p(r))dr (this is the same as L p(r)dr as p £ 1Z); this fixes the 
vertex of the cone of 4> at the origin. 

5.4 Hydrodynamic limit of the delta particle process. 

We now study the hydrodynamic limit for the particle delta process defined in Subsection 15.21 We introduce 
partitions V^ of Z into intervals I^> of length £ where, denoting by L x the interval which contains x, Iq = [0, t— 1] 
(Z?W is now completely specified). We take i equal to the integer part of £ _/3 with /3 <G (0, 1), p <G L°°(M., [0, 1]) and 
(by an abuse of notation) we write 

4P(v) = \ E »?(»), A * ] (p) = ti [ w p{r)dr (5 - 15) 

yetf 5 JeI " 

not making explicit the dependence on e. 

Introduce an accuracy parameter of the form e a , < a < /3; the parameter (3 G (0, 1) is fixed while a will change 
depending on the time. Let 5 £ , a ,^(p)be the set of particle configurations which (e , a, /3)- recognize the macroscopic 
density p £lZ defined by: 

Qs,«Ap) = U G X : Mr,) - l( p )| + ^(77) - n(p)| < e a , sup \A^( V ) - A^(p)\ < e a \, 

L x£Z\{I a Ul L } > 

I = integer part of e~ p , < a < (3 < 1, (5.16) 

where h(rf) and r(t7) are defined in (|2.1j) . and J L is the smallest D*^ measurable interval which contains both h{rf) 
and e~ 1 h(p); I R is defined analogously with reference to 11(77) and R(p). 

Proposition 17. Let p €lZ and t]q ' € Ge,a,p(p), then for any a' £ (0,a) such that 0/ < 1 — /3 and a' < /3/8 the 
following holds: for any k > 1 i/iere are coefficients Ck so that 



P 



e,(5,± /-> / <5,±^ ^ -, „ _fc 



v£7eGe, a >M' ) >l-c fe £ fc (5.17) 



Remark. By iteration the result extends to any finite macroscopic time interval and we can also get a supremum 
over times: 

Corollary 18. Under the same assumptions of Proposition [T7| for any integer m > 1 and for any k > 1 £/iere are 
coefficients cu so that 



[ fl K-™ e S^,^)}] > 1 - <^ fc (5.18) 



n=l 



To show Proposition [17] we need to control the position of the quantiles R a and L 6 of the process evolving with 
the exclusion by the macroscopic time 5e~ 2 . Here a — A e -2 ti b — B € -i t and A and B are Poisson processes of 
parameter ej. These bonds only depend on the exclusion dynamics governed by Lq. 

Sharp convergence of the exclusion process to the solution of the heat equation. Abusing notation denote T t r\ 
the process in X with initial configuration r\ evolving only with the exclusion generator Lq. The SSEP evolution is 
close to a linear diffusion in the following sense: for any n > 2 there is c n so that for any t > 

sup \v(xi,...,x n ;t)\ <c n r n ' s (5.19) 

(:ei,...,:e„)6ZJ 
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where U\ is the set of all n-tuplc of mutually distinct elements of Z, 

n 

v(xi, ...,x n ;t) = E\ ]J[{T t ?7(xj) - u t (x l )} 

and Ut(x) solves the discretizcd heat equation 
dut(x) 1 



Aut(x) = y(ut(x + 1) +ut(x- 1) - 2u t (x)j, mo = V 



(5.20) 



(5.21) 



dt 2 

(|5~T9|) is proved in [2]. The solution of ([531]) is 

V 

Pt(x, y) the transition probability kernel of the symmetric nearest- neighbors random walk. The solution of the heat 
equation starting from p is G t p (recall (|5.10|) V Thus, since n <E G E ,a,p(p) 



\u e - H (x) - G tP (ex)\ < d(et- l ' 2 +e 1 -h- 1 / 2 +e a +et- l ' 2 e-^-^) < c^ 1 "^" 1 / 2 + e a ) 



(5.22) 



The first term takes into account the difference between p t and Gt, the bound follows from the local central limit 
theorem. The second terms takes into account the variations of Gt in an interval of length e 1- " > el. The 
contribution of the difference between averages of 77 and p in good intervals (i.e those not in I h U I R ) is bounded by 
e a ; the contribution of the intervals in I h U I R by rf" 1 / 2 ^^"). 



Boundary control. Using Chcbishcv's inequality and (|5.f9l) 



E 



)aW(T £ - h t)) - A x l) (p t )\ 2n ] < c(e Pn + {e^H- 1 ' 2 + e a ] 2n + [te" 2 ]-"/ 4 ) 
Let 7 > 0, then since 77 £ Ge,a,/3(p), 

J2 P[T e -2 t i 1 (x) = l]<c' k e k , J2 P[Te-HV(x) = 0} < c k e k 



(5.23) 



As a consequence 

By the hypotheses on p, 



x>e~ x — 1 x<— e- 1 — I 

K(T e - HV ) < e" 1 -^; L(T e - H r,) > -e" 1 ^] > 1 - c^ fc 
G t p(r)dr < c' k e k , f dr[l - G t p(r)\ < c' k e 

Jr<-e—* 



I ~k 



(5.24) 



(5.25) 



r>e — ' 



\Wi 



Thus in the event inside the probability in (|5"^|) . for \x\ > e' 1 '^ , \A x '(T e -2 t n) - A x '{G t p)\ < c' k e k . Together 
with (|5.23p this proves that if a' < min{/3, 1 — /3, a, 1/8} and 7 < j3 + a' then for any k > 1 there are coefficients 
Ck so that 



P 



sup 



\A^(T S - HV )- 4P(G tP )\ <e<*' 



> 1 - c k s k 



(5.26) 



Quantile bounds. We now fix t = S and use (|5.26[) to obtain bounds for the quantiles defined in (|5.9[) and (|5.6 
Take 7 < a'/2 and recalling ()5.9j) let 

r' := e- x K j5 {G s p) + £ - l+a '- 2 -< . 



Then 



P(\ Y. T s-*sv(x) - e- x f G s p(r)dr 



<ce a 'e- 1 -~<) > 1- 



Cfc£ 
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On the other hand by the definition of the quantile R^ s (Gsp) and by (J5.25I) 

G sP {r)dr = / G s p{r)dr - / G s p(r)dr - / G sP (r)dr 

eR' Jn^(Gsp) Ju.i s {Gnp) J e~i 

< jS - c'e Q '- 27 
with d = mhi{Gsp(r) : \r — R jS (Gsp)\ < 1} > 0. Hence with probability > 1 — c k e k , 

£ J2 T e- 2 sv{x) < ]5 - de a '- 2 ~ 1 + ce a '£-i < jd - \s a '~ 21 

R'KxKe- 1 -'' 

Let R b (r/) be the quantile defined in (|5.6j) with b = B e -2 S , B t a Poisson process of rate je. Observing that for any 

P(\eB e -2 S -j5\ <s^- K ) >l-c k e k 



we get that for k small enough and with probability > 1 — c k £ k , 

e Y, T E -2 S r,(x) = eB E -* s > jd - s^~ K > 3 5- C -e a '-^ > e £ T e -2 S r,(x) 

x>R b (T E ^2 s v) R'KxKe- 1 -"' 

that implies R b (T E -2 S ri) > R' = £~ 1 R^ s (Ggp) + e~ 1+a ~ 27 . Using an analogous argument for the lower bound we 
get 

P(\eR b (T e -2 S r 1 )-RJ s (G s p)\<e a '- 2 '>) > l-c k e k , 

(5.27) 
P(|eL Q (T £ -2 5?/ ) - \J s {G s p)\ < e a "^) > 1 - c k e k ; 

the second inequality is proved by using the same arguments for a = A e -2$, A t being a Poisson process of rate ej. 

Proof of Proposition [TTl By the definitions (f578|) and (J5TT2|) . 

rf E 'h- = T A ^ B ^T e -2 SV , p 5 s '~ = T^Gsp. 

Since the left and right boundaries after applying T are the quantiles before applying it, inequality f|5. 17[) for the 
delta— processes follows from (|5.26[) and (|5.27|) . The same argument applies for the delta+ processes. □ 

5.5 Hydrodynamic limit of interfaces 

Proof of Theorem 1131 We call r > the time t fixed in Theorem [T3l For each n £ N we let <5 € {t2~"} and 
consider the evolutions ^' ' in a bounded time interval, t <T = 2 N+n 5 = 2 n t, N an arbitrary, fixed non negative 
integer. We have by f)4.13[) . 



= 2eB e -2 (k+1)5 - 2eB e -2 k5 + 2eJ2 % 8 2% +1)S {x) - 2e ]T rf;h%{x) (5.28) 

x>a x>o 

By dUT]) for k = 1, . . . , 2 n+N and x > y, 

eeAtsW - <'-iUv) = 2e(.x -»)-&£ <'->U*) (5-29) 

z=y 

By (|5T8)) and (f5T28|) (f5T29]) ([5^4]) we then get that for any 7 > and any t e {kS : k < 2 N+n } 



lim P 
e-s-0 



sup \ef e 'X (s) - f{{ex)\ > 7 = (5.30) 
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In the next section we shall prove that for any t: 

lim sup|^ 2 ~"' + (r)-^ 2 ~"-(r)| =0 (5.31) 

and that there is a function <p^ (r), r G R, t > 0, so that 

lim sup|0| r) (r) _$*""'- ( r )| = (5.32) 

Then by (lQ0l) - (lQl>(lO2l) 



lim P 

e->0 



sup |e^lL(s) - 4 T \ex)\ > 7 1 = 0, t G T(r) (5.33) 



zGZ 



where 

T(r) = {fc2" T V, fc G N, n G N} (5.34) 

Since r G 7"(r) we have proved (|5.4I) for t = t and since r was arbitrary, Theorem [13] is proved. □ 



6 The macroscopic evolution 

In Subsection 16.11 we prove that as <5 — > the macroscopic delta processes converge -that is, we prove (|5.31|) and 
(|5.32[1 - and that <f> t is well defined by 

61 = lim 6 f ' , t > 

We also collect some properties of the macroscopic evolutions <pt and pt, in particular monotonicity properties of 4>t 
and existence of boundary points for both motions. In subsection 16. 21 we construct macroscopic stationary profiles. 



6.1 Existence and regularity of the macroscopic profiles 

Proof of (|5.3ip and (|5.32[) . Let r > and 8 G {r2 _n ,ri G N}. We shall first prove by induction on fc that for any 
such S, 

Bup|^+(r)-^-(r)|<#. (6.1) 

reK 

(|6.ip holds for fc = because 

o ' + (r) = max{0 o (r), j<5 + \r\}, fi '~(r) = <p Q (r). 

Suppose next that (|6.1[) holds for fc — 1, then by the maximum principle (for the linear heat equation) and calling 
t= (kS)-, 

\$ + (r)-$-(r)\<j6 

hence (|6.ip holds for fc because 

fikt( r ) = maxj^/W, J'( fc + i) 5 + \r\], 4>u ( r ) = max ifiki ( r ), J^ + Ml 

(|6.1|) and (|5.31[) are thus proved. 
It is not difficult to see that 

4>t~ ( r ) — fit ~( r )> fit + ( r ) — fit '"( r ); & = ^' f° r some integer fc > 
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Thus for any n £ N and t > 

^ 2 ""'-(r) < 0[ 2_(n+1) '-(r) < ^ 2 ~"'+(r) (6.2) 

Hence for any fixed t, <f)\ '~ (r) converges pointwisc to a function that we call 4>t v) (it may depend on r). Since 
by definition \<f> t ' (r) — cj> t ' (r')\ < \r — r'\, the convergence is in sup norm and (|5.32p is then proved with c/> t T (r) a 
Lipschitz function with Lipschitz constant 1. □ 

(t) 

In the next Theorem we prove that <jy t is independent of t and also regularity properties of this function. 

Theorem 19. The function <p t (r) in Theorem \13\ is continuous in r andt, more precisely there isc>0 so that for 
all t, t! such that \t — t'\ < 1 and all r and r' , 



\4> t {r)-4>Ar)\<c^\t^¥\, \Mr)~Mr')\<\r-r'\ (6.3) 

Denoting by 8 n any sequence of positive numbers such that <5„+i = <5„/2 then 

cf>t = lim <^ n,± , Vi > 0, (6.4) 

with 4> t n, ~ monotonically increasing and <j> t n monotonically decreasing. 
Proof. Let t > and s > 0, recalling (|5.10[) we have 

G s 4-(r) < ^"(r) < G,$-(r)+J8. (6.5) 

The first inequality is obvious, the second one follows from 

4> S (k+!)s = max {G«50^,j(fc + 1)5 + |r|} < max{G s ^f + j8,j(k + 1)5 + \r\} 
= j6 + G 5 <t> 5 ^, 

because (j>^ > jkS + \r\. 

Since <f> t '~{r) is Lipschitz, it follows from (|6.5p that \4>t+ s ( r ) ~ 4>t~ ( r )l — c^/s + js and, by taking 8 — > 0, 

\<fc\r)-$\r)\ <c^J\t-t'\, for aU r and |* - t'\ < 1 (6.6) 

We shall next prove that ^ r '(r) is independent of r. Obviously <f>j. T (r) = 4>\ T (r) if, recalling (|5.34j) . t' <G T(t) (or 
viceversa). We next suppose that r and r' are not related in such a way. We fix T > and want to prove that 
4>P(r) = 4r'\r). Let 8' = T'2~ n , let 8 < 8' and ft such that kS < 8' < (k + 1)8. Then 

4>8'~ = Toaax{Gs'-kS^ k s~,j8' + \r\} 



< max{G^_ fe5 <^- +j(8' - k5),j5> + \r\} < tff +j(8' - k8) 



because c/)^ s > jk8 + \r\. 

By iteration 4> T '~ < <p^~ + jN8 if N is the cardinality of {k : kS 1 < T}. Thus 

Sl _ £ _ 

and by taking first 8 — > and then 8' — > we get <pf < 4>t ■ The opposite inequality holds as well by interchanging 
(5 and 8' in the previous argument. □ 
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Proof of Theorem [lj From (|5.4[) we get for all 7 > 

IimP(sup|e V Tj <g_ a (x)--{e(b-a)-Meb)-Mea)} > 7) = 



v a<b' 



ea:G[a,6] 



Since </>t is Lipschitz there is pt G -k 1 such that, given any ro G M, 

Mr) = Mro)+ [ (l-2p t (r'))dr' (6.7) 



and since by (|(j.3|) the Lipschitz constant is 1, p t has (almost surely) values in [0, 1] and this proves (|1.2[) . We prove 
in Theorem |2"31 later that p t € 7?., i.e. —00 < L(p t ) < R(pt) < 00. □ 

Proof of Theorem [2} Recall that p\ is the same process as p t '~ . For any t and 5, 

" p 5 t '-(r)dr- I pt{r)dr = <k{b)-<f>\'-{b) + 4'-{a)-<k{a) (6.8) 

J a 

From (J6.4D , for all r and any sequence 5 n so that (5 n +i = J n /2, 

^'"(r)<^(r)<^ 1+ (r) (6.9) 

so that from (J6TT1) and (|6"lfll 

^•"(r)dr - / Pt (r)dr < ^ + (b) - ^-"(6) + t 5 '-(a) - 0*" + (a) < 2j<J n (6.10) 

a 

6.2 Stationary solutions 

We say that a macroscopic interface (p € Vq is stationary if, <po = <p implies (pt = <p + 2jt. A macroscopic density 
p G 7?. is stationary if po = p implies pt = p. Here <^t and pt are the dynamics given by Theorem 1131 and Theorem 
[H respectively. 

If (j> is stationary, then the density p associated to (p via (|6.7p is stationary because by (|6.7|) 

(l-2p t {r'))dr' = (l-2p (r'))dr' 1 for all r and r 

'r J ro 

Viceversa the macroscopic interface <p corresponding to a stationary density p is not necessarily stationary. In fact 
by @2) 

Mr) = Mro) + I (l-2p t (r'))dr' 

Jr 

and since pt = p: 

4>t{r) = (po(r) + [(/>t(ro) ~ <l>o(ro)] 
hence (pt is stationary but only modulo arbitrary time dependent vertical shifts. 
Let 

p(r):=U-2jr, for |r| < £ *(r) := . J + 8 ^' f ''"^ (6.11) 

Theorem 20. The macroscopic interface <p and the associated macroscopic density p are stationary. 

If £q approximates <p i n the sense of f|5.2[) . then the theorem says that the rescaled interface process as seen 
from its vertex converges in the sense of Theorem[13]at any macroscopic time t to the initial value <p shifted by 2jt. 
An analogous statement holds for the particle process (but the stationary density profile does not move). 

Theorem [20] is proven in the next section by introducing a deterministic harness process on M z , a discrete time 
process that approaches <p t and whose stationary solution is directly computable. 
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6.3 Monotonicity 

We collect some monotonicity properties of the macroscopic interface inherited from the microscopic dynamics. We 
tacitly suppose hereafter that the initial data G V , namely that 4>(r) = V (r) = \r\ for all \r\ large enough. 
The following lemma is a direct consequence of the definition of <pt ■ 

Lemma 21. For any t > 0, 4>t(r) > (\r\ + jt), r£l, and 

lim \Mr) - (M + ]t)\ = 

\r\— >-oo 

Proof. Let i = rk2~ n . It follows from the inequality <p T t ~ < 4>t < </>^ that there is i? = i?„ !r so that 

|r| + jt < &(r) < |r| + j(t + 2~ n ), \r\ > R. 

D 



We shall next establish inequalities relating evolutions with different values of j, we thus add a superscript j 
writing <pf' and (j)\ 3 ' ' ' . 

Lemma 22. Let j < j' then for all t > 



4 3) -Jt>4 J ' ] -ft- (6-12) 

Proof. It is enough to prove that 

4>f 5, ~ ] (r) - j<5 > ^"' <5 '" ) (r) - A if <t> > ij 
Let l', r' be such that 

W ' ~ ( r ) = Gsip(r) 7 i/ < r < r', and ipg ' '~ (r) = |r| + j'<5 elsewhere 
By the maximum principle Gg<t> > G^, then 

#'' _) (r) > Gs^r) > $' A -\r), h'<r< r' 

and a fortiori: 

^'-'(r) > 4 J ' A ~Hr) - (/ - j)6, i/ < r < r' 

By definition 0£ 7 ''*' _) (r) > |r| + j<5, relso that 

$' 5 --\r) > \r\+jS = $''*>-\r) (f - j)S, r $ (l',r') 
which concludes the proof. D 

6.4 Existence of boundaries 

Recall the definition in Subsection [573] of the boundaries h(p),R(p) of a density p G 7Z. They are also the boundaries 
of the interface <j> t which corresponds to p t . 

Theorem 23. The boundaries L(p t ), R(pt) of a density p t as defined in Theorerr[l\ starting from p G 1Z are finite. 
In other words, p € 1Z implies pt G H. 

Proof. We shall prove the theorem in the framework of interfaces. We thus want to prove that the boundary points 
of the interface are finite, that is, — oo < L(<% ), &(<$') < oo for initial p G 1Z and all t > 0. Let <f>v > be the 
stationary interface for the j'-evolution. If j' < j is small enough, cf> < (jft > so that, by Lemma |22"1 

4 j) -jt< 4>t 3 ' ] -j't = ij,) . 

This implies that the boundary points of (fr t are bounded by those of <jyj ' and the theorem is proved. D 
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7 The harness process 

We consider now the (deterministic) Harness Process proposed by Hammersley [5] with "moving cone" boundary 
conditions and prove that with the diffusive scaling this process also converges to the macroscopic evolution <j>. 
Let H : Z — >• M. and define the operator by 

(MX .) . *(*-■> + "(*+■> (7 ,, 

but to keep notation light we drop the parentheses and write QH(x). Let (H n (x), x G Z, n G Z + ), H n (x) G M be 
the deterministic process satisfying the discrete heat equation: 

H n+1 (x) := QH n (x) = & n+1 H Q (x) (7.2) 

Here n is time and x is space. 

Duality Let X* be a symmetric nearest neighbors random walk on Z with Xfi = cc and p n (x,y) := P(X% = y) 
the probability that the walk goes from x to y in n steps. Then, a simple recurrence shows that 

H n {x) = Y,Pn(x,y)H (y) = E(H (X%)) (7.3) 

yGZ 

Traveling wave solutions A family of traveling wave solutions of this process are 

H(x) := ax 2 + b, 
where a and b are arbitrary constants. Indeed, 

H n := Q n H = H + 2an. (7.4) 

Moving boundaries For v = (vi, V2) G Z x R, let V v : Z — >• K be the cone defined by K,(aO = \x — V\\ + V2] call 
v the vertex of V v . Let 

%„ := i( : Z 4 I : A'(i) = K,(a;) for all but a finite number of x G Z > 

% := UtjgzxE'Ht) 

For A' G K v define 

l(A) := min{£ G Z : A'(£ + 1) ^ K(^ + 1)} 
r(A) := max{£ G Z : K(l -\) + V v {t - 1)} 

The harness process with moving cone boundary conditions and initial interface Aq G "H(o,o) is defined for n > 1 by 

A^(x) := maxjeAti^)- ^(o,2 m )(x)}. (7.5) 

So that K 3 n G "H (0 ,2.m)- 

A travelling wave solution A of this process is associated to H : if the initial interface is given by 

^t+2jx 2 , \x\<j- 
K&x) := 8J " 4 ; (7.6) 

I \x\, \x\>i 

then the evolution (|7.5[) at time n gives a translation of the initial interface: 

K(x)=K^(x)+2jn. (7.7) 

The sides of the cone y = \x\ are tangent to the parabola y = -^ + 2jx 2 at the points (-jj, ^j) and (jj,jj)- 
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Hydrodynamic limit Let K be a Lipschitz function on W and define 

$^(r):= e (^ 2(] ([ £ -V])) (7.8) 

Proposition 24 (Hydrodynamic Limit). Let <p t be the evolution of Theorem \13\ with initial condition <f> and let 
K^ the evolution (|7.5p with initial interface Kq (x) = <f>(ex). Then, 

lim sup |$j £) (r) - 4> t (r)\ = 0. (7.9) 



t(> 



i el 



where the reseated process $J e (r) is defined in (|7.8p . 

The proof of Theorem (50] follows from the above proposition: 
Proof of Theorem 1201 Taking K e J as the explicit stationary solution in (|7.6p with J = ej and $( £ ) the corre- 



sponding renormalization as in 

(j) t (r) = lim $>\ e \r) = $ e) (0 + 2ji = <j>{r) + 2jt (7.10) 



xi 



where the first identity is consequence of (|7.9p , the second one comes from the stationarity of K given by (|7.7j) and 
the last one is a computation based on the explicit expressions of K^ e > and <fi. □ 

To prove Proposition [24] we introduce the delta processes associated to K . 



The delta harness processes We define the delta harness processes K\' ' , K J e ' as follows. Take S > 1, fix 

(o,o), -^o 15 ^ = max{ii:o' ,5 '~,y (0i( 5 j) ] 



an initial condition /£Tq ' £ "%(o,o)i -^Q) ' = max {^o ' i ^?o,5j)} an( i define iteratively 



K 1 /^ := e^^l/^f, if £e[[nS],[(n+l)S]-l},n>Q 

,5,— r jy-J,5,- 

' xi '■= max A, ;;, 



^Ml : = m^{-^S]_i' V (o,n52 J )}, if n>\ (7.11) 



#, 



n<5] 



= max{if [ ^jl 1 ,V(o,( n+ i)«2j)}, if n > 1. 



Both processes evolve with (j7.ll) in the time intervals [[nS], \(n + 1)6] -ljflZ and update the interface at times [n8]: 
the process delta— takes the max with the cone with vertex (0,ro<52j) while the process dclta+ takes the max with 
the cone with vertex (0, (n + 1)<52j). The following dominating Lemma follows immediately. 

Lemma 25. If K% S '~ < K& < K„' 5 - + , then K 3 /'~ < K\ < K\ A+ , for all £ > 1; for all 6 > 1, J > 0. Furthermore, 

sup {K' i ' s ' + (x)-K 3 /'-(x)}<623. (7.12) 

xez,e>o 

Hydrodynamic limit of the delta processes Take a macroscopic initial condition (j> as in Theorem [T3l Fix 
j > and take J = ej, fix S > and e small such that 5e~ 2 > 1, take 8e~ 2 in the place of 5 and define 

with initial Kq (x) = 4>(ex), which implies 

$^ ± (r,O) = ^(r) = 0( £ [re- 1 ]). 

From Lemma [25l we have 

$^» < *?(r) < $^< + (r) (7.13) 

sup($?'* + (r) -*?•*■>)) <2#. (7.14) 
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Proposition 26 (hydrodynamics). Let 4> be Lipschitz. Then, there exists a constant C > such that, 

supl^'^r)-^^)! < CW-V- 2 *, (7.15) 

r 

for any j5 > 0. 

Proof. Take t < S. In this case <f> t obeys the heat equation. Let W[ be Brownian motion with starting point r; 
then the solution </> t is given by (f>t(r) = G t <fi(r) = Ecp(W[). By duality (|7.3p and assuming X n and W t are defined 
in the same probability space, 

l^f^W-^MI = lE^eXpJ)-^^))! 

< El^eXpJ)-^^)! < E\eX\^-W;\ (7.16) 

< C^e 1 " 2 ^, forany/3>0, t < S, (7.17) 

where in (|7.16|) we used that </> is Lipschitz and in (|7.17|) the dyadic KMT coupling between the Brownian motion 
and the random walk [8], Theorem 7.1. 
At t = 6 we have 

| maxl^l^, % [e -2 5]2j£2) } - max{4^,V {0 ,s23)}\ < C^e 1 " 2 ^ + s, (7.18) 

because the two cones differ at most by e. Changing the constant C, ()7.18|) is bounded by C6^e 1 ~ 2 ^, so that 
iterating (fTTS|l [(< + 1)/S] times we get l|715jl . D 

Proof of Proposition l24l As a consequence of (|7.13[) and (|7.14[) . 

\<Pt-^\r)\ < |^- ^±1 + 1^ -$^1 + 1*^-^1 
< 26 + CtS 1 -^ 1 - 2 ? + 26. 
Taking first e -> and then <5 -> 0, we get (|7T9|) . D 



8 Conclusions 

In this section we summarize the results we have obtained so far. In Theorem 1131 and Theorem [T] we have proved 
convergence in the hydrodynamic limit to a deterministic evolution for the interface and, respectively, the density. 
The limit interface 4>t is Lipschitz continuous in space and continuous in time, see Theorem 1191 At each time 
t > it "belongs" to a cone, in the sense that there are real numbers l(</> 4 ) and R(4>t) so that <j>t{r) = \r\ + jt for 
r ^ (L(4>t) , R.(4>t)) , Theorem l23l The limit particle densities inherit analogous properties from the interface. 

The interface evolution <p t is characterized in terms of a sequence of upper and lower bounds </> t ' which in the 
limit 6 — > become identical. cf> t ' are solutions of time-discrete Stefan problems in the sense that they are obtained 
by alternating linear heat diffusion and motion of the boundaries. 

We miss however a proof that the limit evolution satisfies the Stefan problem described in the introduction for 
the particle density. We do know however that the stationary solution of (|1.7[) - (|1.8p is indeed stationary for the 
limit evolution, Theorem 1201 The formula for the velocity of the boundaries is quite natural once we observe that 
the levels of the solution of the heat equation have velocity —p"/(2p' t ) (p 1 and p" the space derivatives of p). To 
get ()1.8[) we need to add the information that at the endpoints p' = —2j which is consistent with the analysis of 
the stationary solution. 

We have proved in Theorem [4] that there is a stationary measure for the particle process at fixed J = je > 0; we 
miss however a proof that in the limit e — >• it becomes supported by the stationary solution of (|1.7I) - (11.8[) . even 
though this is stationary for the limit evolution. 
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